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Abstract. We investigate the coefficients generated by expressing the falling factorial 
(xy)~ as a linear combination of falling factorial products x-y— for l,m = l,...,k. Al- 
gebraic and combinatoric properties are discussed, some in relation to Stirling numbers. 

1. Introduction 

Let x- denote the falling factorial power, 

x- = x(x — 1) . . .(x — k + 1) , k e N, k > 1 . 

We think of i 6 1, but the definition can make sense in more general settings. Problems 
from discrete Fourier analysis - distant from the topics considered here - led us to falling 
factorial powers of products expressed as 

(i) M k = 
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There is the obvious symmetry c 
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c^j. Since c[l = 1 the interest begins when k > 2. 
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and that the coefficients that do appear are positive. 

Here are the values for cf^ displayed as a symmetric matrix: 
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(k) 

The numbers c\ ) ^ have a number of interesting properties that are the subject of the present 
paper. We found a recurrence relation, several closed-form expressions (which appear rather 
different from each other), a natural combinatorial interpretation in terms of conjoint ranking 
tables, and we can extend these results to products of more than two variables. 

All of this was new to us, but we probably should not have been surprised to discover, 
later, that the combinatorics had essentially been considered in other contexts. To make the 
approach taken here self-contained and more readable we have rederived (briefly) some of 
these earlier results, with references. We also mention several questions that we were unable 
to answer. 

We would like to thank Jiehua Chen, John Gill, Michael Godfrey, and Donald Knuth for 
their comments. 

2. Stirling Numbers and Falling Factorial Powers 

(k) 

One can solve for the coefficients c lr ^ in any particular case, but that there should generally 
be such an expansion emerges from the connection between falling factorial powers, ordinary 
powers and Stirling numbers; see [3], whose notation we follow. 

In combinatorics one defines the (unsigned) Stirling numbers of the first kind by 

k] 

, = the number of permutations of k letters which consist of I disjoint cycles. 



In particular 
(3) 

For us, the important fact is 
(4) 



if k < I. 



x- 



i=i 



x 



Stirling numbers of the second kind, denoted by 



are defined by 



the number of partitions of a set of k elements into I nonempty subsets. 



Here 
(5) 

The special values in 
Corresponding to ([ 

(6) 



1 



if k < I. 



$ and ([S]) will come up in adjusting summation indices, 
one has 



x 



i=i 



x- 



We included the phrase 'Generating Functions' in the title of the paper because the sequences 
\k] fk 
I ' \l 



and ci® each count something and each appears in an expansion in powers, 



([6]), and ([I]), much like classical generating functions. 
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Proceeding with the derivation of (PQ), from and ([6]) we have 

k r , -i k 



p=i l j p=i 



x p y p 
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= Ef- 1 )'-- J E 
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fe p p 

= EEE(-d 1 -' 
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It is now a matter of swapping the summations and using a property of Stirling numbers. 
First, employing Knuth's version of Iverson's convention for sums we write 

k p p 



I [ I m 



x-y- 



EE Ef- 1 )'-" 

p=l i=l m=l 



P I J P I J. 



M m 



l,m,p 



k 
p 



v I I p I „j 



^-[1 < / <p][l <m <p][l <p< Jfe]. 



But now 

[1 < p < k] [1 < I < p] [1 < m < p] — [1 < I, m < k\ [max(Z, m) < p < k] 
and therefore 

fcpp r, -i /■ \ /■ \ k k 



EEEt-D'-" 

p=l i=l m=l 
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Next, since 



for p < I and 



for p < m we can allow the innermost 



Pi ) P V j. 



summation to begin at p — 1 without affecting the result. Thus 

k k k r, -i /• n r n 

i=l m=l p=l 

where 
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i,m=l 
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This establishes the expansion (TjJ and provides a formula for the coefficients. It is not clear 
from this expression that the c's are nonnegative. We will deduce this from a combinatorial 
characterization in Section |H 
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When I = k, 



k 

E 

p=l 
k 

k 



(-1) 



k—p 



k 

m 

Thus the Stirling numbers of the second kind appear in the last row (or column) of the 
matrix of the c's, as we see in ([2]). 

There are a few more arithmetic properties of the cf^ that we wish to note, expressed 

most easily in terms of the symmetric matrix whose (I, m)-entry is c\ ^. Let £W be the 
k x k diagonal matrix whose entries are 

\k] 



J lra 



\k-l 



6, 



lm ■ 



(k) 

Let be the upper triangular matrix with nonzero entries 



OS 



2 Jim 



I < m. 



One can check that 

C« = 4 fe )£( fe )(4 fe) ) T . 

In turn it follows that C^ k > is invertible and has the same signature as Since the 

diagonal entries in T,^ alternate sign, has \k/2] positive eigenvalues and k — \k/2] 



negative eigenvalues. Moreover, the diagonal entries of S% are I's so detS^ 

k 

detC (fc) = JJ(- 1 ) fc ^ 



(k) 



1, and then 



i=i 



For all k it appears that the rows (columns), the diagonals, and the anti-diagonals of 
are all unimodal; this is quite visible for k = 9 in ([2]). Stronger than that, numerical evidence 
suggests all are log-concave. We are only able to show this for the anti-diagonals, 



,(k) 



,(*) Jk) 



• • • ' C i+l,m-l' C l,mi C l-l,m+n ■ ■ ■ ■ 

Recall that a sequence {a^} is log-concave if 

aidi-2 < . 

This is easily seen to be equivalent to the condition 

O'lO'm < a i-i a m+i for all m < I — 2. 
'k 



It is known, [6], that the Stirling numbers 



P 
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are log-concave in / for fixed k. Thus 
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and using this in ([7j) we have 



c 
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which is what we are required to show. 

3. Recurrence 

Falling factorial powers satisfy 

Ax^= (n- l)a£=i, 
where A is the forward difference operator. A little more generally, 

A(arV a ) = (Ax L )y^ + x L {Ay^) + (Ax l -)(Ay^) 
= lx l -^-y^ + mx l -y— + lmx l -=±y^ 

Then, on the one hand, 

fc 

A{xy)~ = c< i%( lx —y- + ™h— + lmx l -^y r -^) 



l,m=l 
fc-1 fc 



k fc-1 



fc-1 



/=0 m=l i=l m=0 

On the other hand, apply the identity 



l,m=0 



v- = -(v— + kv- 

V 



with v = (x + l)(y + 1) to write 

A(xy) i =((x + l)( 2 /+l))^-(x 2 /) 
1 



((* + l)(y+ l))^ + k((x + \){y + - {xyf 



and then 



(x + l)(y + l) 

A(Tv] k _ V" ( fc +i) (x + l^y + l)- A ( fe) {x + lfjy + l)^ 



r ( fc ) r L,2i 
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Comparing the two expressions for A(xy)~, and rearranging terms, we have shown: 
Theorem 1. The coefficients cf^ satisfy the recurrence 

cffiSUi = 4% + (l + + (™ + l)c£ +1 + ((I + l)(m + 1) - k)cf^ m+l . 

We comment that other proofs of the recurrence are possible, for example one that uses 
the recurrence relations for Stirling numbers. 

4. Combinatorial Characterizations 

Conjoint analysis is a method in marketing that assigns weightings to independent at- 
tributes of a product. It is only the first step that we consider, that of setting up a conjoint 
ranking table. Suppose there are / > 1 choices for one attribute (e.g. price) and m > 1 
for the other (e.g. color). Form an I x m table, in which each cell represents the two 
attributes considered jointly; hence the contraction 'conjoint.' Let k < Im and rank the 
conjoint preferences from 1 to k. Since we allow k < Im not every pair of attributes need 
have a ranking, but we do insist that every individual attribute must enter into the ranking 
at least once. That is, every row or column must have at least one ranked cell. In particular 
k > max{/,m}. One can consider the remaining cells as left blank or filled in with zeros. 

For a given I and m we refer to such an object as a k- conjoint ranking table. This is 
quite a general concept and it is easy to imagine many examples. Here is one other. A 
graduate student is signing up to take qualifying exams, and is required to fill out an I x m 
table indicating preferences. Each of the I rows refers to a subject area, such as algebra, 
combinatorics, analysis, etc.. Each of the m columns refers to a particular professor; we 
assume that they are all equally capable of asking about any of subjects. The student will 
be asked k questions, where k < Im, and may put the numbers 1 through k in any of the 
cells, indicating preferences for who asks what kinds of questions. However, the student 
is not allowed to avoid any professors, and must put at least one number in each column. 
Similarly, the student is not allowed to avoid any subject areas, and must put at least one 
number in each row. 

The recurrence relation allows us to derive a combinatorial interpretation and characteri- 
zation of the c's as counting the number of such tables. 

Theorem 2. The number of k-conjoint ranking tables of size I x m is l\m\ c|^. 

Before giving the proof we observe the following consequence. 
Corollary 1. The satisfy 

c l m = f or Im < k and cf^ > for Im > k . 

This is so because when Im < k, the conjoint ranking table is too small to fit all k numbers 
and so cj^ = 0. However, when Im > k, there are enough cells in the table to be ranked 
from 1 to k. Furthermore, there must exist at least one way of placing the numbers so as 
to satisfy the row and column constraints since the conditions 1 < l,m < k imply that 
k > max(/,m). Thus cf^ > 0. 

One can see the pattern of zeros for k = 9 in ([2]). The parabolic shape of the boundary 
between the zero and nonzero coefficients becomes more pronounced as k increases. 
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Proof of Theorem [H To prove the theorem we first define an equivalence relation on the 
set of fc-conjoint ranking tables, namely two tables of the same size are equivalent if one is 

(k) 

obtained from the other by permuting the rows, or columns, or both. Let d\ ^ be the number 

of equivalence classes. We show that af^ = cf^ by showing that the same recurrence relation 
obtains. 

The following four cases are properties of any table in an equivalence class: 

(1) A; is alone in both its row and column; 

(2) k is alone in its row but not its column; 

(3) k is alone in its column but not its row; 

(4) k is neither alone in its row nor its column. 

These cases are illustrated in Figure [TJ We count the number of equivalence classes by 
computing how much each case contribute to d\ 



(k) 



(k-1, 1-1, m-1) 



Case (1): 
"k" is alone in 
both its row and 
its column 



(k-1, 1, m-1) 



Case (3): 
"k" is alone 
in its column but 
not its row 



(k-l,l-l,m) 



Case (2): 
"k" is alone in 
in its row but 
not its column 



(k, 1, m) 



Case (4): 
"k" is neither 
alone in its row 
nor its column 



Figure 1 . Diagrams illustrating the four cases in the proof of Theorem [2j 

Case (1): If we remove both the row and column containing k, then we are left with 
a (k — l)-conjoint ranking table of size (I — 1) X (m — 1). Hence this case contributes 
d (k ' 1] 

u l-l,m-V 

Case (2): If we remove the row containing k then what remains is a (k — l)-conjoint 
ranking table of size (I — 1) x m. Since the number k could have been in any of the 
m slots of the row just removed, the contribution is md^ 1 ^. 

Case (3): The analysis is identical to that of Case (2) and contributes ld\ k ~^ v 
Case (4): Whereas in previous cases we removed slices of a fc-conjoint ranking table, 
now we will replace a zeroed entry of a (k — l)-conjoint ranking table of size I x m 
with k. There are lm — {k — 1) zeroed entries in such a table, so this case contributes 
ilm-{k-l))df- l) . 
Adding the contributions from each of the four cases yields 

dtl = dtiXi + mdl'% + ld^l\ + (Im - (k - l))d^ 
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which is the same recurrence as for cf^ (reindexed) . Considering small tables verifies that 

(k) (k) 

d]^ starts out like c]^, and the two quantities are therefore equal. 

Since we allowed for permuting I rows and m columns, the number of A;-conjoint ranking 
tables of size Ixmis l!m! cf^. 

□ 

Here we make contact with earlier work, for the core of the proof above is counting a set 
of binary matrices, and these have been counted in other ways; see Chapter 9 in |2j. One 
approach is to use the inclusion-exclusion principle. We would like to give this argument to 
show how it leads to another expression for cf^ (which also appears in older literature, but 
not as interpreted here). 

Fix / and m, take k < Im, and let O be the set of binary matrices of size I x m with l's 
in exactly k positions. Then 

Now let k > max{Z, m} and let C be the subset of O which have at least one 1 in every row 
and in every column. We want to find \C\. 

For the index i running from 1 to I let Ai be the subset of O whose z'th row has all O's, 
and for the index i running from 1 to m let Ai+i be the subset of O whose i'th column has 
all O's. Then 

C = ^n^n-n A l+m , 

where 

Ai O \ Ai • 

By the inclusion-exclusion principle 

l+m 

\c\ = \o\-J2 lAiHAJ- IAnA 2 nA 3 |+' • -+(-i) i+m |An^2n- ■ -rvl* +m | . 

i=i h<h n<«2<*3 

To compute the general sum 

^l^n-.-nAJ 

we have to have to distinguish the matrices that have zeroed rows from those that have 
zeroed columns. Suppose among the h sets A^ , . . . , Ai h that p of them have zeroed rows. 
Then h — p have zeroed columns. For a fixed p there are ( J ways to select p rows to zero out 
and there are (^.) = ( m _™ +p ) ways to select h — p columns to zero out. After these choices 
there are pm+(h—p)l—p(h—p) zeros total, and so there remain lm—pm — (h—p)l+p(h—p) = 
(I — p)(m — h + p) cells amongst which we place k ones. There are then 

'(/ — p)(m — h + pY 
k 



ways of doing this. Thus 



p=0 



l\ { m \ ( (I — p)(m — h + p) 
p J \m — h + p) \ k 
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and 

^ * ,_^ h {l\( m \ f(l-p)(m-h+p) 



' ' \p J \m — h + p J \ k 

h=0 p=0 KF/ v 1 ' v 

Multiplying \C\ by k\ distinguishes the nonzero elements, whether they are k distinguished 
balls tossed into bins or the numbers from 1 to k. Then dividing by l\m\ allows for permuting 
the rows and columns. The end result is evidently the same as counting the number of k- 
conjoint ranking tables, and hence 

(g) c ( fc ) _ kl \n\ - kl V^V^/ i\hf l \f m \({l-p){m-h + p) 



^ l ' m /!m! |W| l\m\ ^ ^ x ^' \p)\m — h + pj\ k 

There is one more approach and one more formula. The equation (OQ) can also be written 
in the form 

(9) ( x v) = Y b W(*\(V 



where 



k J ^ l ' m \l 

l,m=l 



-« _ l] - m [ c (k) 



m 



l,m y,\ l,m ' 

We understand the binomial coefficient to be defined for nonintegral x by 

'x\ T(x + 1) 



kj k\T(x - k + 1) ' 
and we use 

r(x + 1) 

x{x-l)...(x-k + l) - 



r(x- k + 1) 

(k) 

We see that the b\ ^ give a count of the number of the I x m binary matrices with exactly 
k ones such that each row and column has at least one 1. We have a further comment on the 
combinatorics of Q) when x and y are integers. The left-hand side, n?) is just the number 
of ways to select k elements from anixt/ array. How does this jibe with the right-hand 
side? 

For any given selection of k cells there is a unique minimal subarray (smallest number 
of rows and columns) such that every row and column of that subarray contains a selected 
cell. This is illustrated in Figure [2J On the left, the darkened cells correspond to a selection 
of k = 8 cells, and the arrows indicate the rows and columns of the subarray in which the 
selected cells are contained. On the right we see the subarray extracted, and note that every 
row and column contains a selected square. Counting all possible such subarrays thus counts 
the ways of choosing k elements from the big array. This is what the right-hand side in 
equation (jUJ) does, for rfj is the number of ways to choose anlxm subarray, and then 
for each such subarray we count the number of ways to populate it with k entries such that 
no row or column is void. That multiplier is precisely of^. 
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column 
indices 



V V V V 



row 
indices 



1 2 3 4 5 6 



1 

2 
3 
4 
5 
6 
7 
8 
9 



2 3 5 6 



2 
3 
4 
5 
7 
8 



selection of 8 entries 
from 9x6 array 



corresponding subarray 
in which every row and column 
contains selections 



Figure 2. Combinatorial Interpretation of equation 
From Theorem [TJ the 6's satisfy the recurrence 



(10) 



(1 + !)("»+ + + Ci,J + ((* + !)(«» + 1) - 

But explicitly in P it was noted that Mobius inversion applied to Q implies 

in 6« = ]T(-i; 



■,(*) 



(,(*) 



(,(*) 



Lm 



,t=l 



Furthermore, a similar formula seems also to be found in [I]. Mobius inversion is kin 
matrix inversion and it is interesting to see how the latter can be used to derive 
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For k 2 pairs j/j 



unknowns b\ k ^. In matrix form, A/3 = £ where 



k 2 , we treat Qj as a system of linear equations for the k 2 



T) ft 1 ) 

?) to 



(a) ft 



r 2 2 ) ft 



2 



fi') 



I 2 / 



(T)0r) Mfif) ••• (T)(T) (T)( y f) C¥)CW 



xi\ yi 

kJ\k 

X2\ (V2 

k)\k 



Of) (T) 



.4 



b {k) 
1,1 

6 (fc) 
"1,2 



6 (fc) 
U l,k 

°2,1 

6 (*) 
"2,2 



u(k) 

°L{i),M(i) 



J k,k 



k ) 

k J 



k ) 



L ( Xk T) J 



Here L(z) and M(i) are mappings that convert from linear indexing to matrix indexing: 
i - 1 



L(i) 

To go the other way 
Thus, succinctly, 

.4 



k 



+ 1, M(z) = ((i-l) modJfe) + l, i = l,...,A; 2 . 
i = k(l — l) + m, l,m = 1, . . . , k . 



k 



M \L(i)J\M(j)J ' W * 

We want to choose the to make A invertible, and solve for A -1 . We find that if we 
take = L(i) and = M(i), so that 

L(i)\/M(i)\ 



^ \L(j)J \M(j)J 
then A is invertible with 

(12) det A = 1 , Ar 1 = (-l) L ^+ L U)+M(i)+M(j) ^ 
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This closely mimics the well-known phenomenon exhibited by the Pascal matrix and its 
inverse: 

In [5] we prove a more general result so we will not give the details behind (j!2p here. Briefly, 
the first step is to observe that A is the element-wise product (Hadamard product) of the 
two matrices A^ 1 and A^ y \ where 

iJ W,)J- ^ \M(j)J> 

and A^ is, as a whole, block lower-triangular with k x k size blocks while in A^ y ' the separate 
k x k blocks are each lower triangular. See Figured It follows that A is lower-triangular. 





A {1 



A {y 



A 



Figure 3. Shapes of A^ x \ A^ and A. The matrix A on the right is the 
Hadamard product of the matrices on the left. 

The diagonal elements of A are all 1 so det A — I. Defining B = (-l) L »+ L 0')+ A/ «+ A/ (i) A id 
one then shows directly that BA is the identity matrix, and this requires some special 
identities for binomial coefficients. 

From the formula for A~ l let us derive ( ITT]) . We have 



P = A- 1 £. 



Lm 



Pi , where i — k(l — 1) + m, 



or 



■,(*) 



i=i 
k 2 



k ) 



3=1 
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For L(i) and M(i), using 1 < m < k it is easy to verify that 



L(i) 



i - 1 



+ 1 



k(l - 1) +m- 1 



1 = Z 



M(i) = ((i-1) mod A;) + 1 = ((k(l - 1) + m - 1) mod fc) + 1 = m. 
To simplify the L(j) and M(j) terms, express j as the combination 

j = ks + t, s G {0, l,...,fc-l}, t e {1,2,..., A;}. 
Then just as above, 



L{j)=L(ks + t) 



ks + t — 1 



s + 



t - 1 



k 



+l=s+l 



M(j) = M(ks + t) = ((ks + t - 1) mod jfe) + 1 = (t - 1) mod k + l = t- l + l = t. 
Hence, 

fc 2 



i=i 
fc-i fc 



y^^_^L(i)+L(fcs+t)+M(i)+M(fcs+t) 

EB- 1 



L(i) \( M(z) 



L(ks + t)J \M(ks + t) 

s + l)t 
k 



L(j)M{ks + 1) 
k 



\l+s+m+t 



s=l i=l 

which is (TTTT). 



Z\ / m\ / st 



We conclude this part of the discussion by noting that we have three different expressions 
for the c's (or the fe's), ((7|), (jHi) and (TTTT) - and we do not know algebraically how to derive 
one from another! 

5. More than two Variables 

For an expansion 

(xxX 2 ■ ■ ■ X n )~ = ^ ■ ■ ■ Xn— , L = (Zi, Z 2 , ■ ■ ■ , l n ) , 1 < k < k , 

L 

of a product of more than two variables all the results in the preceding sections have natural 
extensions. Modifications to the earlier arguments are straightforward and so we record the 
outcomes with little additional detail - the chief problem is notation. We follow the generally 
accepted conventions on multi-indexing. In particular 

L! = Zx!Z 2 ! ■■■l n \. 
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The formula for the c's in terms of Stirling numbers is 



c (*) - \ \ i 

or simply 



p=i l -i i=i 



p 



4" = E(-D'-' fil {£ 

p=i L^J L 



if we allow ourselves the analog to the the multi-indexed case of binomial coefficients and 
write 

'l}-fl{l 

' i=i 

There is a natural extension of the product rule for the forward difference operator and it 
can be applied just as before to obtain a recurrence relation. It will pay to invest in a little 
extra notation. We write 

L + 1 = (Z f + 1: % = 1, . . ,,n) 
and for a subset S C {1, . . . , n} we write 

Ls = (h'- i € S) , L s + l = (h + l:ieS), Lg = (l,:ie{l,...,n}\5). 
The result is 

n 

C i+1 ) = C L ^ _ ^ C L+1 + £ £ ^11 ^) C is+l,% ' 

m=l |5|=m ie5 

The generalization of a 2-dimensional conjoint ranking table allows for n independent 
attributes (color, price, shape, . . . ) with U choices for the z'th attribute. Based on the 
recurrence one can then show that Llc^ is the number of ways to fill in an l\ x I2 x • • • x l n 
conjoint ranking table with the numbers 1 through k, insisting, as before, that each attribute 
must enter into the ranking at least once. Again this implies the nonnegativity of the c's, 
more precisely 

cf = if Y[h < k and cf > if Y[h > k. 

i i 

Extensions of the alternate formulas (|S]) and (ITT]) are more complicated to write. For the 
former, to keep the final result from being too cluttered we use the notation 

and for a multi-index M = (jni,m2, ■ ■ ■ , m n ) with 1 < rrii < h and ||M|| = h we let 

n n / 

m) = n u - £ 1 i r H(i s - m s ) n i s 

i=l r=l \ s<r s>r 

Then 



c 



/i=0 M,||M||=h v 7 j'=l v J 



While the derivation of this formula is only an extension of the argument for two variables, 
some discussion will help make the form clearer. Let "slices" refer to the higher-dimensional 
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generalization of the notion of rows (or columns) for matrices. As before, we use inclusion- 
exclusion to count , the number of distinct l\ x . . . x l n (0, l)-tables with exactly k ones and 
no zeroed slices; that's the formula without the factorials in front. In the outer summation, 
the index h is the number of zeroed slices. The inner summation then runs over possible 
ways to distribute the h zeroed slices across the n different dimensions, rrii is the number of 
zeroed slices in dimension i and is the number of ways to select a particular set of rrii 
slices to zero out. After counting ways to zero the slices, we must then count the number 
of ways to fill the remaining cells in the table with k ones. The total number of cells in the 
array is Y\d=i h an d *&(L, M) counts the remaining cells by subtracting the number of zeroed 
slices from the total number of cells. Then (*^^ ) is the number of ways to distribute k 
ones among these remaining cells. 

The analog of the formula pip based on matrix inversion is developed in [!5j . It reads 
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